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Abstract 
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magnitude of the NLO QCD corrections is about 10% in most of the parameter space and can 
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2j +'rjet+ l^r signature from the and the decays including the NLO QCD effects, and find 
an observable signal at a 5(j level in some parameter region of the minimal supergravity model. 
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I. INTRODUCTION 

The Higgs mechanism [l| plays a key role for the understanding of the spontaneous elec- 
troweak symmetry breaking in both the Stand Model (SM) and the Minimal Supersymmetric 
Stand Model (MSSM) Searching for Higgs bosons is one of the most important missions 
for the upcoming CERN Large Hadron Collider (LHC). The MSSM contains five physical 
Higgs bosons: two neutral CP-even bosons and H^, one neutral CP-odd boson A^, and 
the charged boson pair. The is the hghtest and SM-like Higgs boson, while the others 
are non-SM-like ones whose discovery will give the direct evidence of new physics beyond 
the SM, especially charged Higgs boson. 

At hadron colliders, the charged Higgs bosons if^ could appear as the decay product of 
primarily produced top quarks if the mass of is smaller than nit — nib- For heavier H^, 
single charged Higgs boson production associated with heavy quark, such as gh — > H^t j^, 
qh — > q'bH~ [4], and qq, gg thH^ are the main channels for single charged Higgs boson 
production. The channels for pair production are qq annihilation and the loop-induced gg 
fusion process jol. These processes have large production rates, but also suffer from large 
QCD backgrounds, especially when the mass is larger than mt + mf,. Another attractive 

n 

channel is single charged Higgs boson production associated with W boson [7] . The dominant 
partonic subprocesses at the LHC are bb W^H^ at the tree-level and gg —>■ W^H"^ at the 
one-loop level Sj. For the bb annihilation process, the supersymmetric electroweak (SUSY- 
EW), the 0{as) pure QCD and the supersymmetric QCD (SUSY-QCD) corrections have 
been calculated in Ref. 9| lOj 111, respectively. In this paper, we use the dimensional 
reduction (DRED) [l^ scheme to regularize both the ultraviolet (UV) and the infrared 
(IR) divergences while in Ref. lOj] the gluon was given a finite small mass to regularize IR 
divergences. We will focus on the case of /i > which is favored by the recent measurement 



of the anomalous magnetic moment of the muon 



13|, /i is the Higgs superfield mass term in 



the superpotential, so the SUSY-QCD corrections are relatively small and can be neglected 
as shown in Ref. . For simplicity, in our calculations, we neglect the bottom quark mass 
except in the Yukawa couplings. Such approximations are valid in all diagrams, in which the 
bottom quarks appear as initial state partons, according to the simplified Aivazis-Collins- 
Olness-Tung (ACOT) scheme [l4j. Moreover, we only consider the process bb —>■ H^W^ 
since the cross section for the process hb — H^W~ is the same if we choose all the relevant 
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parameters to be rea . 

n 

Recently, in Ref. [15| the authors investigated the viabihty of observing charged Higgs 
bosons produced in association with W bosons at the LHC at LO level, using the leptonic 
decay H —>■ r^Vj. and hadronic W decay. In this paper we also give the Monte Carlo 
simulation results of the above signal, but in the minimal supergravity (mSUGRA) [l^ 
scenario including the NLO QCD effects. 

The arrangement of this paper is as follow. In Sec. II, we show the LO explicit expressions. 
In Sec. Ill, we present the details of the calculations for both the virtual and real QCD 
corrections. In Sec. IV, we give some analysis on the signal and background. Sec. V are the 
numerical results for total and differential cross sections and the Monte Carlo simulation 
results. Sec. VI contains a brief conclusion. The relevant coupling constants and the lengthy 
analytic expressions are summarized in Appendix. 



II. LEADING ORDER CALCULATIONS 



The tree-level Feynman diagrams for the subprocess b{pi)b{p2) — > H {p3)W~^{pi) are 
shown in FiglH and its LO amplitude in n = 4 — 2e dimension is 



= Mj^^ + M^o + M^^, + M^^ , (1) 



with 



= " 2M /it mn ^^'''' - '"''"(ZMs + M3)], (2) 

where Sa = sin a, Cq, = cos a, S/3_q, = sin(/? — a), c^_q, = cos(/3 — a), = tan/?. Mandelstam 
variables s, t, and u are defined as follows: 

s = ipi+ P2)^ t = {pi- psf, u = {p2- Pzf. (3) 

Mj's are reduced standard matrix elements, which are defined by 

Mi{2) =v{p2) Ap^)Pr{l)u{pi), 
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^3(4) =V{P2) ^4 /iPi)PRiL)U{pi), 

^5(6) = v{.P2)Pr(l)u{pi)pi ■ e{pi), 
^7(8) = v{p2) APr(l)u{pi)pi ■ e(p4), 

^9(10) = viP2)PR{L)U{pi)p2 ■ e(p4), 

and 

Mii(i2) = v{p2) APr{l)u{pi)p2 ■ e(P4), (4) 

with the projectors Pl,r = (1 =F 75)/2 • 

The LO total cross section at the LHC is obtained by convoluting the partonic cross 
section with the parton distribution functions (PDFs) G^g/p in the proton: 

cj-^ = y dxidx2[Gb/p{xi, fXf)Gi/p{x2, fif) + (xi ^ X2)]a^ , (5) 

where /i/ is the factorization scale and = J ^^|M'^p(ir is the Born level cross section 
for b{pi)b{p2) —* H~{p3)W~^{p4), in which the colors and spins of the outgoing particles have 
been summed, and the colors and spins of the incoming ones have been averaged over. 

III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO QCD contributions to the associated production of H~ and through bb 
annihilation process consist of the virtual corrections, generated by loop diagrams of col- 
ored particles, and the real corrections with the radiation of a real gluon or a massless 
(anti)bottom quark. For both virtual and real corrections, we use DRED scheme to regu- 
larize all the divergences. 

A. Virtual corrections 

The virtual corrections to bb H~W~^ arise from the Feynman diagrams shown in 
Figl2l which consist of vertex, self-energy and box diagrams. We carried out the calculation 
in 't Hooft-Feynman gauge and used the dimensional reduction in n = 4 — 2e dimensions to 
regularize the ultraviolet, soft and coUinear divergences in the virtual loop corrections. In 
order to remove the UV divergences, we use the modified minimal subtraction (MS) scheme 
to renormalize the bottom quark mass and wave function, while for the top quark mass 
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and wave function we use both the MS scheme and the on-shell (OS) scheme and compare 
them. Denoting nibo, ^to, V'w and ipto as the bare quark masses and the bare wave functions, 
respectively, the relevant renormalization constants 6mb, 6mt, 5Zi,l,r and 5ZtL,R are then 
defined as 

rriM = mb + 6mb, 
mo = mt + Smt, 

iPbO = (1 + SZbL)'^^^bL + (1 + SZbR)'^^^bR, 

i^to = (1 + SZtLf'^^tL + (1 + 6Z,RY'^^tR. (6) 

with ip{b,t)L = PL4'(b,t) and ip(b,t)R = PR4'(b,t)- After calculating the self-energy diagrams in 
FiglJl we obtain the explicit expressions for all the renormalization constants as follows: 



— — — — dL-i^A, 

TUb rrit At: 

rrit Air \ m"^ / \e 

6Z2ln)t = -^Cf [Re(So + B,) - 2m?Re(i?^ - B[)] {ml 0, m?), (7) 

where A = 7^ + ln(47r), Cp = 4/3, and B' = dB/dp^, Bi are the scalar two-point 

integrals [l7|. 



The renormalized virtual amplitude can be written as 

= M""""" + M^""". (8) 

Here M""^*^" contains the radiative corrections from the one-loop vertex, self-energy and box 
diagrams, as shown in Figl2], and M™" is the corresponding counterterm. Moreover, M™'"'^"' 
can be separated into two parts: 

/ 

j^unren ^ M"" + M^°^' , (9) 

a=a 

where a denotes the corresponding diagram indexes in FigJ2J Using the standard matrix 



5 



elements from Eq.(jl]) they can be further expressed as 



12 



1=1 

12 

^box = J2fl'°''Mi, 



(10) 



1=1 



where and fi°^ are the form factors, which are given exphcitly in Appendix. The coun- 
terterm contribution M^""" is separated into M^°^^'^^ and M^°"*^*\ i.e. the counterterms for s 
and t channels, respectively, which are given by 



j^con 

j^con{s) 

J^con{t) 



rrih 



Ml 



con{t) 



m 



2\2 



{ [mt{t - ml)6Zt - {ml + t)5mt] 



m: 



conit) 



aM2 - b[{t - m])5Zt - 2mt5mt]{M^ + 2M9)} 



with 



errit 



erribtp 



(11) 



(12) 



y/2MwSwt/3 ' \/2MwSw 
After adding all the terms above, the renormalized amplitude is UV finite, but still 

contains the IR divergences, and is given by: 



IR 



27rr(l - 2e) 



4 + ^ ) 

e 



(13) 



with 



-Cf, A'l=~Cp. (14) 
Here the IR divergences include both the soft and the collinear divergences. The soft diver- 



gences are canceled after adding the real emission corrections, and the remaining collinear 
divergences can be absorbed into the redefinition of PDF [18], which will be discussed in the 
following subsections. 



B. Real gluon emission 

The Feynman diagrams for the real gluon emission process b{pi)b{p2) — > {p2,)W^ (pi) + 
g{p5) are shown in FigJS] 



The phase space integration for the real gluon emission will produce both soft and 
coUinear infrared singularities, which can be conveniently isolated by slicing the phase space 
into different regions defined by suitable cutoff parameters. In this paper, we use the two- 
cutoff phase space slicing method , which introduces two small cutoffs to decompose the 
three-body phase space into three regions. 

First, the phase space can be separated into two regions by an arbitrary small cutoff 6s, 
according to whether the energy {E^) of the emitted gluon is soft, i.e. < 6s\^/2, or 
hard, i.e. E^ > 6sy/s/2. Correspondingly, the partonic real cross section can be written as 

a^ = a^ + a^, (15) 

where and cr^ are the contributions from the soft and hard regions, respectively, 
contains all the soft divergences. Second, in order to isolate the remaining collinear diver- 
gences from , we should introduce another arbitrary small cutoff, called collinear cutoff 
6c, to further split the hard gluon phase space into two regions, according to whether the 
Mandelstam variables satisfy the collinear condition —6cS < ui^2 = (pi,2 — V^Y < or not. 
Thus, we have 

where the hard collinear part a^*" contains the collinear divergences, while the hard non- 
collinear part a^'^ is finite and can be numerically computed using standard Monte-Carlo 
integration techniques and can be written as 

da^ =-Y}MydV^. (17) 

Here (iFs is the hard non-collinear region of the three-body phase space. 

In the next two subsections, we will discuss in detail the soft and hard collinear gluon 
emission. 

1. Soft gluon emission 

In the soft limit, i.e. when the energy of the emitted gluon is small, with E^^ < 6s\^/2, 
the matrix element squared ^jM-^p for the process b{pi)b{p2) — > {p3)W^ (p^) g^p^) can 
be simply factorized into the Born matrix element squared times an eikonal factor $ejfc- 

^|M«(66 ^ H-W+ + g)f (47ra,/i,2^)^|M^|2$e.., (18) 



7 



where the eikonal factor $ejfc is given by 



Moreover, the phase space in the soft hmit can also be factorized as 



(19) 



dT^{hh + g) dT2{bb H-W+)dS. 

where dS is the integration over the phase space of the soft gluon, which is given by [1 



(20) 



dS 



\3-2e 



dE^El '^^dVL2~2e- 



2(27r) ^0 

Hence, the parton level cross section in the soft region can be expressed as 



(21) 



(22) 

Using the approach of Ref. jl9|], after analytically integrating over the soft gluon phase space, 
Eq. fl22l) becomes 



with 



r(l-6) [An^4X^ (A^ ^^1^4. 
27rr(l-2e) V s ) \ \e^ ^ e ^ 

= 2Cf, Al = -ACf In 5,, Al = ACf In^ 5, . 



(23) 



(24) 



2. Hard collinear gluon emission 

In the hard collinear region, i.e. > 6sy/s/2 and —6cS < -Ui 2 < 0, the emitted hard 
gluon is collinear to one of the incoming partons. As a consequence of the factorization 



theorems 



20[ | , the squared matrix element for bb ^ H + g can be factorized into the 



product of the Born squared matrix element and the Altarelli-Parisi splitting function for 



(bb) b{b)g [2JJ 
Y\M^{bb 



22|,i.e. 



9. 



12 collinear 



B\2 



-2Pbb(z) -2Pr,(z 



+ 



bb\ 



(25) 



ZUi ZU2 

where z denotes the fraction of incoming parton 6(6) 's momentum carried by parton b{b) 



with the emitted g. 



uon taking a fraction (1 — z), and Pij{z) are the usual Altarelli-Parisi 



splitting kernels 2l|]. Explicitly, 

P,t,{z) = Piiiz) 



CF^^ + CF^il-z). 

1 — z 2 



(26) 
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Moreover, the three-body phase space can also be factorized in the coUinear hmit, and, for 
example, in the limit —6cS < ui < it has the following form 

dTsibb + g) dr2{bb R-W"-- s = zs)—^^ -dzdui[{z - l)ui]-\ 



IQtt^TH - e) 



(27) 



Here the two-body phase space should be evaluated at the squared parton-parton energy zs. 
Thus, the three-body cross section in the hard collinear region is given by 13] 



da 



as r(i-6) 

27rr(l - 2e) 



2\ 



1 



5/ [Pbb{z)Gb/pixi/ z)Gi/p{x2) 



+Pll{z)Gi/p{xi/z)Gb/pix2) + (xi ^ X2)] — 



dz f 1 



dxidx2, 



(28) 



where G^iypi^x) is the bare PDF. 



C. Massless (anti) quark emission 

In addition to the real gluon emission, a second set of real emission corrections to the in- 
clusive production rate of pp H^W^ at the NLO involves the processes with an additional 
massless (anti) quark in the final states: 

gb bH~W+, gb bH-W+. 

The relevant Feynman diagrams for massless (anti) quark emission (the diagrams for the 
antiquark emission are similar and omitted here) are shown in FigjH 

Since the contributions from the real massless (anti) quark emission contain the initial 



3 



state collinear singularities, we also need to use the two cutoff phase space slicing method 
to isolate those collinear divergences. Because there is no soft divergence in the splitting of 
g — > bb, we only need to separate the phase space into two regions: the collinear region and 
the hard noncollinear region. Thus, according to the approach shown in Ref . [l^ , the cross 
section for the processes with an additional massless (anti) quark in the final states can be 
expressed as 

^^add ^ J2 a^{a(3 ^ + X)[G^/p{xi)Gp/p{x2) + (xi ^ X2)]dxidx2 

{a=g,(3=b;b) 



as r(l - e) / Anjj,; 



2 \ n 



-^)^c' [Pb9{z)Gg/p{Xi/ z)Gi/p{x2] 



277 r(l - 2e) 

dz (\ - z\ 

+Gb/p{xi)Pig{z)Gg/p{x2/ z) + (xi ^ X2)] — i 1 dxidx2, (29) 
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where 

P,,{z) = n,{z) = ^-[z' + {1 - zn (30) 

The first term in Eg. (129 p represents the noncollinear cross sections for the two processes, 
which can be written in the form: 

rfa^=^^|M"^|W3, (31) 

where a and /3 denote the incoming partons in the partonic processes, and dV^ is the three- 
body phase space in the noncolhnear region. The second term in Eg. (1291) represents the 
colhnear singular cross sections. 

Moreover, the top momentum in FigJl](c) and (e) (as well as in the corresponding b 
emission Feynman diagrams) can approach the top mass shell, which will lead to a singularity 
arising from the top propagator. Following the analysis shown in Ref. j23|, this problem 
can easily be solved by introducing the non-zero top width Fj and regularizing in this 
way the higher-order amplitudes. However, these on-shell top contributions are already 
accounted for by the LO level tH~ and iW~^ productions with a subseguent decay, and thus 
should not be considered as a genuine high-order correction to H~W~^ associated production. 
Therefore, to avoid double counting, these pole contributions will be subtracted in our 



numerical calculations below in the same way as shown in Appendix B of Ref. 



23|. 



D. Mass factorization 

As mentioned above, after adding the renormalized virtual corrections and the real cor- 
rections, the partonic cross sections still contain the collinear divergences, which can be 



absorbed into the redefinition of the PDF at NLO, in general called mass factorization |18 |. 
This procedure in practice means that first we convolute the partonic cross section with the 
bare PDF Ga/p{x), and then rewrite Ga/p{x) in terms of the renormalized PDF Gn/p(x, Uf) 



In the MS scheme and DRED scheme, the scale dependent PDF Ga/p{x, iJ,f) is given by [19| 

^ / N ^ / \ v^/ Ixfas r(l - e) /47r/i^ VI dz ^, , , , , 

Ga,/p[x,Hf) =G^/p{x) + 2^(--) _ 1^^;;^ j \ J —PapKz)Gp/p[x/z) 



9^E/ ^PL,{x/y)G^„{x), (32) 

Jx y 
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where P\ are the regulated sphtting functions and PlAz) are the usual Altarelli-Parisi 



0/3 

splitting kernels 



21|, exphcitly 



Cf 
1 



1 + 



PI^[z) = -Cf{1-z)+Cf^5{1-z) 
Pi^{z) = -z{l-z). 



(33) 



After replacing the bare PDF by the renormalized MS PDF and integrating out the coUinear 
region of the phase space defined in the two-cutoff phase space slicing method [lol, the 
resulting sum of Eg. (1291) and the collinear part (the second term) of Eq. ( |28l) yield the 
remaining 0{as) collinear contribution as: 



coll 



a 



as r(l-6) / 47r/i^ ^ ' 
27rr(l-2e)V s 



a=b.b 



Af{a ag) 



+ Al\a ag) 



Gb/pixi, fif)Gi/p{x2, iif) 



a. 



+ {xi ^ X2)}dxidx2 — —^Cf<J^, 

2tt 



(34) 



where 



with 



(6 ^ bg) = Al'ib hg) = GF{2\nSs + 3/2), 



/l- = /irin(-), 



—Gp/p{x/y,fif)Pap{y) 



(35) 
(36) 

(37) 



Papiy) = PaM in(4^4) - PLpiy)- 

y fj'f 



(38) 



The NLO total cross section for pp H W'^ in the MS factorization scheme is obtained 
by summing up the Born, virtual, soft, collinear and hard noncollinear contributions. In 
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terms of the above notations, we have 



cr 



NLO 



dxidx2{ 



Gb/p{xi, fif)Gi/p{x2, fif) + (xi ^ 0:2) 



a + a + a + a )} + cr 



coll 



+ / dxidx2 



Ga/p{Xl, fif)Gp/p{x2, fif) + {Xi ^ X2) 



(39) 



We note that the above expression contains no sing ularities, for 2AX + = and 2AY + 
Al + Al%b bg) + Al%b bg) = 0. Namely, all the l/e^ and 1/e terms cancel in a^^'^. The 
apparent logarithmic 6s and 6c dependent terms also cancel with the the hard noncollinear 
cross section a^'-' after numerically integrating over its relevant phase space volume. 



IV. MONTE CARLO SIMULATIONS 



Based on the work of Ref. jl5|], we discuss the same signal in the mSUGRA scenario 
including the NLO QCD effects. In the signal channel, H~ decays leptonically, H r^Ur 
and decays hadronically, — qq'{q = u, c, q' = d,s). For simphcity, we only consider 
hadronic decays of the r lepton, r ^ Ut- + hadrons. The resulting signature is 2j + rjct+ i>T-, 
where the missing transverse momentum "^j^ is carried away by the two neutrinos, Tj^t comes 
from r decay and 2j come from W boson decay. The transverse mass is defined as: 



m± = \ 2pT,Tjet - cos(A0)], (40) 



where A0 is the azimuthal angle between pT,Tjet and As pointed out in Ref. [l5l], the 
m± distribution will have a peak with the upper edge of the peak given by the mass of the 
charged Higgs boson. The two light jets can be distinguished by calling them hard (with 
momentum phj) and soft (with momentum psj) according to the larger and smaller value 
of their transverse momentum p^, respectively. Our study is performed at parton level, 
without considering parton showering or hadronization, and the detector effects also not be 
considered. Event generation is performed with help of PYTHIA v6.206 [24] and TAUOLA 



v2.7 



25 



261] is used to perform the decay of r lepton. 
The cuts we have used are shown in Table [H which are the same as in Ref. [l5| in order to 
compare our results with theirs. Here the basic cuts define a signal region that corresponds 
to the sensitive region of a real detector and the additional cuts are used to suppress both 
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Basic cuts 


Additional cuts [all in GeV] 


\VTiJ < 2-5 


PT,r,et > 50, > 50 


\7]A < 2.5 

1 IJ \ 


70 < m,-,- < 90 

J J 


ARjj > 0.4 


m± > 100 


ARr^^,, > 0.5 


PT,hj > 50,pT,sj > 25 


PT,jet > 20GeV 





TABLE I: Basic cuts and additional cuts used 



background and detector misidentifications. The dominant irreducible SM background for 
our signature 2j + Tjet+ /Pt comes from W + 2j production with W tv^-. We use 



ALPGEN 27| to repeat the background calculations of Ref . [ISj] , and the same results can 
be obtained. The + 2j background mainly comes from ug and dg initial states, while 
W~+2j background is mainly due to dg and ug initial states. Detailed descriptions about the 
backgrounds and cuts can be found in Ref. |l5| , and our simulation results will be discussed 
below. 



V. NUMERICAL RESULTS 

The arrangement of this part is as follow. First, we present the NLO QCD calculations 
of both total cross sections and differential cross sections. Then we turn to the simulation 
results under several groups of cuts and mSUGRA parameters. 

A. NLO cross section calculations 



In the numerical calculations, we used the following set of SM parameters 28| : 



aewimw) = -- = 1/128, mvF = 80.40GeV, = 91.1876GeV, 
47r 

mt = 174.2GeV, mf,(mf,) = 4.2GeV, as{Mz) = 0.1176, = 0.23122. (41) 



The running QCD coupling as{Q) is evaluated at the two-loop order [29| and the CTEQ6M 



PDF [301] is used throughout this paper to calculate various cross sections, either at the 
LO or the NLO. As for the factorization and renormalization scales, we always choose 
f^f = "^av = {^H± + ^w)/'^ aiid fir = fif, unless specified otherwise. Moreover, as to the 
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Yukawa couplings of the bottom quark and top quark, we took the running masses mb{Q) 
and mt{Q) evaluated by the NLO formula |3l| : 



rribiQ) = U(i{Q,mt)U5{mt,mb)mb{mb), 

miQ) = UG{Q,mt)mt{mt), (42) 



with mb{mb) = 4.2GeV 28|]. The evolution factor Uf is 



Uf{Q2,Q 



d(f) r 



^ ^ asiQi) - as{Q2) jif) 
Air 



dU) = 12 ^ 8982 - 504/ + 40/^ 

33-2/' 3(33 -2/)2 ' ^^""^ 

where / is the number of the active hght quarks. We use both the MS and the OS renor- 
malization scheme for top quark in our calculations and find good agreement in these two 
schemes. We will only show the numerically results in the MS scheme unless specified 
otherwise. 

The values of the MSSM parameters taken in our numerical calculations were constrained 
within mSUGRA, in which there are only five free input parameters at the grand unification 
(GUT) scale. They are mi/2, rriQ, Aq, tan/5, and the sign of /i, where mi/2, rriQ, Aq, /i are, 
respectively, the universal gaugino mass, scalar mass, the trilinear soft breaking parameter, 
and the Higgs superfield mass term in the superpotential. Given those parameters, all the 
MSSM parameters at the weak scale are determined in the mSUGRA scenario by using the 



program package SPHENO 321 ]. 



In Fig. O we show the dependence of the NLO QCD predictions on the two arbitrary 
theoretical cutoff scales 6s and 6c, introduced in the two-cutoff phase space slicing method, 
where we have set 6c = 6s/100 to simplify the study. The NLO total cross section can be 
separated into two classes of contributions. One is the 2-^2 rate contributed by the Born 
level, and the 0{as) virtual, soft and hard collinear real emission corrections, denoted as 
(T^, (T^, a^, and cr^°'' in Eq. fl39|) . Another is the 2 — > 3 rate contributed by the 0{as) hard 
noncoUinear real emission corrections, denoted as o"^*^ and o"*^ in Eq. fl39p . As noted in the 
previous section, the 2 — 2 and 2^3 rates depend individually on 6s and 6c, but their 
sum should not depend on any of the theoretical cutoff scales. This is clearly illustrated in 
Fig. El where <Jnlo is almost unchanged for 6s between 10"'* and 10~^, and is about 25.6fb. 
Therefore, we take 6s = 10~^ and 6c = 5^/100 in the numerical calculations below. 
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Fig. [6] shows the total cross sections for pp H~W~^ at the LHC in both the MS scheme 
and the OS scheme as a function of rriH- for tan/9 = 5,20 and 40, respectively, assuming 
mo = 150GeV, and Aq = 300GeV. The results in the two schemes are almost the same. The 
total cross sections decrease with the increasing rriH- ■ In general, the NLO QCD corrections 
enhance the total cross sections for small tan j3, but reduce for large tan j3. 

In Fig. [71 the total cross sections for pp H~W^ at the LHC are plotted as a function 
of tan/9 for two representative values of mi/2- When tan/5 ranges between 5 and 45, tuh- 
varies from 290 GeV to 185 GeV, and from 595 GeV to 402 GeV for mi/2 = 160 GeV and 
400 GeV, respectively. From Fig. [7] we can clearly see that the total cross sections increase 
with the increasing tan/5 and decrease with the increasing mi/2- For large tan/5(> 40) and 
mi/2 = 160 GeV, the LO and the NLO total cross sections can be over 100 fb. 

Fig. [S] gives the dependence of the K factor (defined as the ratio of the NLO total cross 
sections to the LO ones in the MS scheme) on mn- for H~W^ production, based on the 
results in Fig. [61 It can be seen that the results in the two schemes are in good agreement. 
For instance, the difference of the K factors in the two schemes is within 4% for tan (3 = 5 and 
less than 2% for tan/5 = 20 and 40. In general, the K factor decreases with the increasing 
mn-- For tan/5 = 5, the K factors can increase to 1.1 when rriH- < 400 GeV. While for 
tan/5 = 20 and 40, the K factors decrease below 0.9 when rrijj- > 500 GeV. 

Figs. [HI shows the dependence of the total cross sections for pp —>■ H^W^ production at 
the LHC on the renormalization scale (/i^) and the factorization scale (/U/), with /x^ = /i/- 
We defined R as the ratio of the cross sections (LO, NLO) to their values at central scale, 
fir = fi'f = ^av = {^H± + "^vk)/2, always assuming fij. = fi'f for simplicity. For three values 
of tan (3, the scale dependence of the NLO total cross sections reduced when going from LO 
to NLO in both the MS and the OS scheme. For example, in the MS scheme, the ratio R 
at the LO vary from 0.78 to 1.03 when fir = fif ranges between 0.2mav and 5mav, while the 
NLO ones vary from 0.98 to 1.08, for tan/5 = 40. 

Fig. Uni shows the differential cross sections as the functions of the transverse momentum 
Px of the H~ and the in the associated production of the H~W^ pairs at the LHC. In 
case (1), the NLO QCD corrections can enhance and reduce the differential cross sections 
in the medium pt region of the and the respectively, and are negligible small in 
both the high and the low px region. In case (2), the NLO QCD corrections reduce the 
differential cross section significantly in the medium p^ region of the H~ , otherwise the 
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NLO QCD corrections can be neglected. 

In Fig. [11] we display the differential cross sections as the functions of the invariant mass 
Mh-w+ of the H~W^ pairs produced at the LHC. In case (1), the NLO QCD corrections 
enhance the LO differential cross sections more, which can reach 10%, in the medium values 
of the MH-y/+, but are negligible small in both the high and the low values of the Mh-w+- 
In case (2), the NLO QCD corrections reduce the LO differential cross sections in the high 
values of the MH-y/+, while the corrections are relatively small in the low values of the 
Mh-w+- 

Note that our numerical results of the NLO QCD corrections to the total cross sections 
are different from the ones given in Ref . |l^ , where the corrections are always negative and 
the magnitude can reach 30%. We also used the same parameters as in Ref. [lOj] to compare 
with their results, but our results are still different from theirs. 



B. Simulation results 



Our simulation results for the relevant distributions are shown in Figs. [T2lfT6| which 
include the distributions of the my, the pt for all jets and the missing transverse momentum 
for the signal and backgrounds after the basic cuts, assuming: (1) ttlq = 200GeV, mi/2 = 
147 GeV,Ao = 200 GeV,tan/3 = 50 and fi > 0; (2) mo = 320 GeV,mi/2 = 400 GeV, = 
300 GeV, tan/3 = 50 and fi> 0. In case (1) of those figures, the NLO QCD corrections can 
be neglected for all the distributions, but in case (2), the NLO QCD corrections reduce the 
LO results significantly, which can reach above 10% in some region of the distributions. It 
can be seen that the additional cuts introduced at the LO still work well when including the 
NLO QCD effects. 

In the following calculations of the total cross sections the additional cuts are used. 
Moreover, an integrated luminosity of 300 fb~^ and a r detection efficiency of 30% are taken 
to calculate the significance S/ \/~B. The total cross sections for the backgrounds from the 
final state + 2j and W~ + 2j are about 32 fb and 25 fb, respectively. Now, we add 
the cross sections of the H~W^ and the H~^W~ production together, as well as for the 
backgrounds. Tables [III and IIIII show some representative results of the cross sections and 
the significance, where we can see that the significance can reach above 20 for tan/3 = 50 
and rriH- = 175 GeV. 
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tan/3 = 50. 







Integrated 


cross section 


(fb) 




Parameter 




Signal 


Background 




s/Vb 


mjj- = 175GeV 


LO 


17.6 


57 




22.1 




NLO 


17.2 


57 




21.6 


niH- = 345GeV 


LO 


2.12 


57 




2.66 




NLO 


1.84 


57 




2.31 


niH- = 630GeV 


LO 


0.34 


57 




0.43 




NLO 


0.28 


57 




0.35 


Its under several groups 


of nifj- value, assumin; 


y rriQ 


= ^0 






Integrated 


cross section 


(fb) 




Parameter 




Signal 


Background 




S/^/B 


tan (3 = 20 


LO 


0.46 


57 




0.58 




NLO 


0.42 


57 




0.53 


tan 13 = 40 


LO 


4.60 


57 




5.78 




NLO 


4.34 


57 




5.45 


tan /3 = 50 


LO 


17.7 


57 




22.3 




NLO 


16.7 


57 




21.0 



200 GeV, and 



TABLE III: Results under several groups of tan (3 value, assuming ttiq = Aq = 200 GeV, and 
mi/2 = 150 GeV. The three tan/3 values correspond to mH± = 282 GeV, 224 GeV and 178 GeV 
from top to bottom, respectively. 



Figs. [T7] and [18] show the dependence of the cross sections on m^^- mass and tan/3, 
respectively. In general, the NLO QCD corrections reduce the cross sections in most of the 
parameter space, and their magnitude can be larger than 10%. The horizontal lines in the 
figures correspond to the total cross sections required for S/ ^/B = 5. It can be seen that at 
the LO the signal can be detected at a 5a level for tan/5 > 40 and mi/2 = 150 GeV, and for 
100 GeV < rriH- < 250 GeV and tan /3 = 50, respectively. And the NLO QCD corrections 
have small effects on the above results. Some of our results are different from those given in 



Ref. 



151] mainly due to the difference between the mSUGRA scenario and the one used in 
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Ref. 
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VI. CONCLUSIONS 

In conclusion, we have calculated the NLO QCD corrections to the inclusive total cross 
sections of the H'^W^ pairs produced at the LHC through bb annihilation in the MSSM. 
The NLO QCD corrections can either enhance or reduce the total cross sections, but they 
generally efficiently reduce the dependence of the total cross sections on the renormaliza- 
tion/factorization scale. The magnitude of the NLO QCD corrections is about 10% in most 
of the parameter space and can reach 15% in some parameter region. Finally, we give some 
discussion on the H'^W'^ 2j + Tjet+ signal including the NLO QCD effects, and find 
an observable signal at a 5a level in some region of the mSUGRA parameter space. 
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APPENDIX 

In this appendix, we give the relevant Feynman rules and the form factors for the virtual 
amplitude. First we give the relevant Feynman rules. 

1. h^{H^) -b-b: Aii2)mi, 

^ 2cf3MwSw' ^ 2ci3MivSw' 
where a is the mixing angle in the CP even neutral Higgs boson sector. Here we use the 
abbreviations Sq, = sin a, = sin (3 and so on. 

2. A'^ - b-b: Asnib'j^ 

2MwSw 
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3. h\H^,A^)-W~ -H+ : ^i(2,3)(p. 



I — — , -J 2 



hO{HO,AO) - Ph+, 



2sw 2sw 2sw 

Here we define the ingoing four-momenta to be positive. 
4.b-i~H+: giPL + Q2PR 

iedijUit ieSij-rUbtp 



\/2MwSwtp'' ''^ \/2MwSw 
Here and below, we assume tlie third generation CKM matrix element Vtb equal to 1. 
5. h-t-W- : Gs-f^'PL 



-le 



\plsw 



Below we collect the explicit expressions of the nonzero form factors in Eq. ( fTOl) . For simplic- 
ity, we introduce the following abbreviations for the Passarino-Veltman two-point integrals 
-Bj(,), three-point integrals Ciu\ and four-point integrals -Dj(j), which are defined similar to 



Ref. [171] except that we take internal masses squared as arguments: 







(s,0,0). 




= 


(t,0,m2). 




= Ci{j 


(M^_,0,t,m?,0,0), 




= Ci{j 


(M2„0,t,m2,0,0), 




= Cii^j 


(0,3,0,0,0,0), 




= Ci(^j 


(0,t,M|_,0,0,m2), 


fie 


= Cii^j 


(0,t,M2„0,0,m2), 


A(i) 


= Di^j 


(M2„s,0,t,M^_,0,m2, 0,0,0) 



Most of the above functions contain IR singularities. Since all the Passarino-Veltman inte- 
grals can be written as a combination of the scalar functions Aq^Bq^ Cq and Dq^ we present 
here the explicit expressions for the Cq and Dq functions, which contain the IR divergences 



and were used in our calculations 

:2 3 

a 



Co 



t 

-^ln^(mi - M]j-) + ^In^im^ -t) + lnsln^^^ 



Ml. 



In ^ + Li 



t 



Ml- 



Ml. 



m: 



-Li 



mi 



Ml- 
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Co 



s{t - ml) 
- In s In - — 



Ml- 



w 



(t-m2)2 



[ml ~ Ml^){ml - M^) 



mt 



M^_)(m: 



-Li 



1 + 



Ml_){m\-M, 



smt 



-ln(mi - Ml-) In 



mr 



Ml. 



2Li ( 1 + 
1 

~ 2 



mr 



mt 



2Li 1 + 



m: 



Ml- 



mt 



- In^ ( 1 + 2 In s In 



mt 



ln(m 



mt 



M^)ln 



t — ml 



mt 



2 

mt 



where = (47r/i^/s)T(l — e)/r(l — 2e). For diagrams(a)-(f) in FigJH we get the form 
factors as following, respectively, 



£a ra £a j?a R 

h-h-h- ho - ^^^^i^cpMwsi^{s 



m 



fb fb fb f 

Jb ~ Je ~ Jd ~ J: 



cb 
'10 



a. 



e^mhCaSf3_a 



4:71 cpMwslyis 



m 



HO, 



— ~f& — ~fio — CpSij 
a. e 



a. 



e^mbtp 



27r V2sw{t - ml)^ 



4:7iMwsly{s - m\o) 
{2mta){tB\-mlB^), 



[{C^, + C'i + C^,)s-2B^] 
[{C^, + Cl + C^,)s-2B^] 
[(Co^ + C^ + C2^)s-2i?, 



OJ ) 



/a = V' = 

/2 = Cp^ij 



a. 



a.. 



27r ^/2sw{t 



m 



2\2 



-2h)[ml{B\-?,B',)+t{B', + B\] 



47r ^swit 



mt 



a. 



47r ^swit 



{amt) [45o' + 2{Ml- - 1)0^ + 2(2M^- - t)Ci] , 

-26) [4^0^ + 2{Ml- - t)Cf + 2(2M|_ - ml)C^] 



mt 



^'-^""^'''^TT^s^it-ml) 



fi 



fi 



f'^-^^^^U^r^s^it-ml)^''' 
where a,b are abbreviations for 



a. 



47r ^/2sw{t 

a. e 



mt 



47r ^/2sw{t - ml 
a. e 



(2am,) [2Co^o + ^^(^2 + ^^22) + (M^ " t)iC1 + C^^)] 
(-26) [2Co^o + M^{C^, + CI,) + (M^ - t){Cl + €{, 
{-Ah) [2Cl, + Ml{Cl + C^, + CI, + CI,) - mlCl] 



a 



emt 



V2MwSwtp ' 
For the box diagram(g) in Figj^l we find 



embtf3 
V2MwSw ' 



ft = C^%^-^(-2m,a) { [-2D00 - M^{Du + /^is + + D3 



-Ml- {D2 + Ds + D22 + D23)] - {Di2 + Drs + D23 + ^33)^ - (^^12 + Dq + + D2)u} , 



box 



CpSij 



47r ^/2sw 

£box K ^ 



box 



Cp^ij 



a. 



47r 



{-2b)Dss, 

(46) [-CS + M^{Do + Di + D2 + D3) - ^^2^] , 
(-2mta) [-2(Do + Di + 2D2 + + + D22 + 1^23)] , 
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ft = '^^^'^^^/^^^^^ - Co" + M^{Do + D,) + {t + u- Ml.)D^ + + u)D^] , 
fll^ = CF6,~-j^{-2mta) [-2(^2 + + D,^ + i^ia + ^^22 + 2^3 + i^ss)] • 
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FIG. 1: Leading order Feynman diagrams for bb ^ H . 
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FIG. 2: One- loop virtual Feynman diagrams for bb ^ H 
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FIG. 3: Feynman diagrams for the real gluon emission contributions. 
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FIG. 5: Dependence of the NLO total cross sections for the H production at the LHC 

on the theoretical cutoff scale 6s with 5c = 5^/100, assuming the mSUGRA model with niQ = 
150 GeV,mi/2 = 300 GeV,^o = 300 GeV,tan^ = 40 and ^ > 0. 
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FIG. 6: Dependence of the total cross sections for the H production at the LHC on niH- , 
assuming mo = 150 GeV, and Ao = 300 GeV. 
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FIG. 7: Dependence of the total cross sections for the H production at the LHC on tan/3, 
assuming mo = 150 GeV, and Aq = 300 GeV. 
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FIG. 8: Dependence of the K factor for the H production at the LHC on m^-, assuming 
mo = 150 GeV, and = 300 GeV. 
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FIG. 9: Dependence of the ratio R on the renormalization/factorization scale (//^ = l^f) for the 
H^W^ production at the LHC, assuming: (a) mo = 150 GeV, 771,5^/2 = 170 GeV,Ao = 300 GeV, 
and tan/3 = 5; (b) mo = 150 GeV,mi/2 = 300 GeV,Ao = 300 GeV, and tan/3 = 20; (c) mo = 
150 GeV, mi/2 = 300 GeV, Aq = 300 GeV, and tan /3 = 40. 
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FIG. 10: Differential cross sections in the transverse momentum (pt) of the H~ and the bosons, 
for the H~W^ production at the LHC, assuming: (1) mo = 200 GeV,mi/2 = 180 GeV, = 
250 GeV, and tan/3 = 40; (2) mo = 150 GeV,mi/2 = 400 GeV, = 300 GeV, and tan/3 = 40. 
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FIG. 11: Differential cross sections in the invariant mass {Mh-w+), for the H production 
at the LHC, assuming: (1) mo = 200 GeV,mi/2 = 180 GeV,^o = 250 GeV, and tan/3 = 40; (2) 
mo = 150 GeV, mi/2 = 400 GeV, ^o = 300 GeV, and tan (3 = 40. 
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FIG. 12: niT distributions for the signal and backgrounds after basic cuts, assuming: (1) tjiq = 
200 GeV,mi/2 = 147 GeV, and Ao = 200 GeV; (2) mo = 320 GeV,mi/2 = 400 GeV, and 
Aq = SOOGeV. The left axis scale is for the cross section of the signal, while the right one is for 
the backgrounds. 
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FIG. 13: Transverse momentum distributions of the r jet for the signal and backgrounds after 
basic cuts, assuming: (1) mo = 200GeV,mi/2 = 147 GeV, and Aq = 200 GeV; (2) mo = 
320 GeV, mi/2 = 400 GeV, and Aq = 300 GeV. The left axis scale is for the cross section of the 
signal, while the right one is for the backgrounds. 
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FIG. 14: Missing transverse momentum distributions for the signal and backgrounds after ba- 
sic cuts, assuming: (1) mo = 200 GeV,mi/2 = 147 GeV, and = 200 GeV; (2) mo = 
320 GeV, mi/2 = 400 GeV, and Aq = 300 GeV. The left axis scale is for the cross section of 
the signal, while the right one is for the backgrounds. 
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FIG. 15: Transverse momentum distributions of the hard light jet for the signal and backgrounds 
after basic cuts, assuming: (1) mo = 200 GeV, mi/2 = 147 GeV, and = 200 GeV; (2) mo = 
320 GeV,mi/2 = 400 GeV, and = 300 GeV. The left axis scale is for the cross section of the 
signal, while the right one is for the backgrounds. 
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FIG. 16: Transverse momentum distributions of the soft light jet for the signal and backgrounds 
after basic cuts, assuming: (1) mo = 200 GeV, mi/2 = 147 GeV, and Aq = 200 GeV; (2) mo = 
320 GeV, mi/2 = 400 GeV, and = 300 GeV. The left axis scale is for the cross section of the 
signal, while the right one is for the backgrounds. 
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FIG. 17: H"^ mass dependence of the integrated cross section after all cuts, assuming mo = = 
200 GeV. 
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FIG. 18: tan/? dependence of the integrated cross section after all cuts, assuming mo 
200 GeV. 
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